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Abstract 

The instability of an infinite swept attachment line boundary layer is 
considered in the linear regime. The basic three-dimensional flow is shown to 
be susceptible to travelling wave disturbances which propagate along the 
attachment line. The effect of suction on the instability is discussed and 
the results suggest that the attachment line boundary layer on a swept wing 
can be significantly stabilized by extremely small amounts of suction. The 
results obtained are in excellent agreement with the available experimental 
observations. 
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INTRODUCTION 


The laminar boundary layer which forms on a long cylinder whose axis is 
inclined relative to the oncoming flow is known to exhibit instability to 

small amplitude disturbances at large Reynolds numbers. Depending upon the 
flow conditions these disturbances may take the form of Tollmien-Schlichting 
waves, stationary cross-flow vortices or, if the surface has regions of 
concave curvature, Taylor-Gortler vortices. In general, as the Reynolds 

number is increased, the instability occurs first in the regions of adverse 
pressure gradient on the lee side of the cylinder. Further increases in 

Reynolds number cause the location of the instability to move progressively 

forward towards the attachment line. It follows that the Reynolds number at 
which the attachment line flow exhibits an instability to small amplitude 
disturbances represents a limit beyond which the flow over the cylinder is 
unstable everywhere. 

The motivation for studying the flow over inclined cylinders lies in the 
fact that the results find direct application in situations of engineering 
significance, notably to the flow over high aspect ratio swept back wings. In 
recent years the attempted development of laminar flow wings has resulted in 
renewed interest in the instability mechanisms which can affect a three- 
dimensional laminar boundary layer. For conventional wings operating at 
typical long range cruise conditions a laminar boundary layer would be 
unstable over a large proportion of the surface area and, consequently, wholly 
laminar flow cannot be maintained without the use of some form of artificial 
stability augmentation. The most convenient form of stabilization is that 
provided by distributed surface suction. Surface suction affects the boundary 
layer in two ways. Firstly, the viscous layer is thinned with a corresponding 
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reduction in local Reynolds number. Secondly, and more importantly, the 
vorticity distribution within the boundary layer is modified in such a way 
that a more stable flow is established. The combination of these separate 
effects makes surface suction a particularly efficient stabilization mechanism 
with transpiration velocities of order .1% of the velocity at the edge of the 
boundary layer being sufficient to maintain aerofoil boundary layers in the 
laminar state. Therefore, an investigation of the stability characteristics 
of the attachment line boundary layer with and without surface transpiration 
constitutes a problem of considerable practical importance. 

The significance of the attachment line with respect to laminar flow 
control has been recognized for many years (see Pfenninger (1977)). The 
previous work on this problem has been almost wholly experimental though, 
following Gregory, Stuart, and Walker (1955), there have been many theoretical 
investigations of crossflow instabilities which are important outside the 
attachement region. The small amplitude stability problem was first 
considered by Gaster (1967) who introduced controlled acoustic disturbances 
into the attachment line flow formed upon a faired circular cylinder model. 
From his measurements he concluded that the swept attachment line boundary 
layer was stable to small amplitude disturbances for momentum thickness 
Reynolds numbers up to at least 170. In a later investigation Cumpsty and 
Head (1969) reported that no instability was observed for Rq up to 240. The 
first experimental investigation to detect amplified attachment line 
disturbances was conducted by Pfenninger and Bacon (1969). Their test 
configuration consisted of a 45° swept back wing with a blunt-nosed aerofoil 
section which also had provision for surface suction through closely spaced 
slots. The results of this experiment will be discussed in depth later in 
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this paper. More recently, an investigation of the stability characteristics 
of the attachment line flow on a swept cylinder has been carried out by Poll 
(1979 and 1980). In these tests, which were similar in some respects to those 
conducted by Pfenninger and Bacon, amplified disturbances were observed and 
recorded. By combining the results of all the available experimental data 
including the little known work of Carlson (1966), Poll argued that the upper 

limit for the stability of the laminar swept attachment line boundary layer 

flow occurs at a momentum thickness Reynolds number of approximately 230. 

In this paper we investigate the instability of a three-dimensional 
boundary layer which is directly relevant to the experimental configurations 

of Pfenninger, Bacon and Poll. The relevance of the flow will be immediately 

obvious when we describe it but we postpone until the last section of this 
paper a more detailed discussion of the relationship between the experimental 
situation and our idealized model. 

We consider the flow adjacent to the infinite flat plate defined by 
y " 0 with respect to Cartesian coordinates (x,y,z) with the x and z 
axes lying in the plate. The velocity components in the x and z 

directions are zero at the plate and approach the values U Q j and W Q 

respectively when y + 00 , The normal velocity component takes on the 

constant value Vq at the wall and grows linearly with y when y + ®, The 

important simplifying feature of this flow is that it corresponds to an exact 
solution of the Navier-Stokes equations so that it is not necessary for us to 
make the boundary layer approximation when deriving the basic flow. Thus it 
is sensible for us to discuss the stability of the flow at finite Reynolds 
number and the concept of a critical Reynolds number is tenable. This is in 
contrast to the situation with Blasius flow where a self-consistent 
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perturbation investigation of the linear stability problem does not lead to a 
critical Reynolds number. The latter problem has been discussed by Smith 
(1979) and Bodonyi and Smith (1981). 

The basic flow described above is susceptible to centrifugal instabilities 
because of the curvature of the streamlines in the x-y plane. In particular 
when Wq = 0 the plane stagnation point stability equations of Hammerlin 
(1955) are recovered when a disturbance periodic in z is imposed on the 
flow. Hammerlin showed that a continuous spectrum of neutrally stable 
wavenumbers exists for these equations. The eigenfunctions corresponding to 
these modes decay algebraically with y but it is yet to be shown that they 
are physically relevant in a more realistic situation where they must be 
matched onto a disturbance outside the attachment region. There also exists a 
discrete spectrum of damped eigenvalues having eigenfunctions which decay 
exponentially when y -*• » and it is this type of eigenvalue on which we will 
concentrate our attention here. 

If Wq is nonzero the spanwise velocity component is susceptible to 
Tollmien-Schlichting waves but we show that there is no rational approximation 
which leads to the Orr-Sommerfeld eigenvalue problem. Instead we find that 
the appropriate eigenrelation is of sixth order in y and reduces to that 
obtained by Hammerlin in the limit Wq 0. 

The numerical solution of the eigenvalue problem is nontrivial because of 
the rapidly varying nature of the eigenfunctions and two independent schemes 
were used in our calculations. The most efficient scheme was a fourth-order 
compact finite difference method based on that described by Malik, Chuang and 
Hussaini (1982). In addition the Riccati method was also used to provide a 
check on the eigenvalues obtained from the finite difference method. 
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Our calculations showed that at a critical value of Wg the attachment 
line boundary layer is unstable to discrete Tollmien-Schlichting 
instabilities. Moreover, these modes are connected to the discrete modes of 
the plane stagnation point stability problem found by Wilson and Gladwell 
(1979) so that we have an additional check on the numerical work. The 
eigenfunctions associated with the discrete spectrum change from being of 
typical Taylor-Gortler type at small values of W Q to being of Tollmien- 
Schlichting type when instability occurs. 

The results which we obtain for the basic flow described previously can be 
applied to the experiments of Pfenninger and Bacon (1969) and Poll 
(1979,1980). We shall compare our predicted frequencies with those given by 
these authors and there seems to be little doubt that the instability 
mechanism which we investigate is responsible for that found experimentally. 

The procedure adopted in the rest of this paper is as follows: in Section 
2 we derive the basic flow and formulate the eigenvalue problem governing the 
stability of this flow to spanwise periodic disturbances. In Section 3 we 
discuss the boundary conditions at infinity which must be imposed before the 
eigenvalue problem can be solved. In Section 4 we discuss the numerical work 
required to calculate the basic flow and solve the eigenvalue problem. In 
Section 5 we discuss the asymptotic limit of large suction whilst in Section 6 
we discuss the relevance of our results to experimental observations. 


2. THE BASIC FLOW AND THE DISTURBANCE EQUATION 

We consider the flow of a viscous incompressible fluid of kinematic 
viscosity v adjacent to a flat plate defined by y = 0 with respect to the 


- 6 - 


Cartesian coordinate system (x,y,z). The y axis is taken to be normal to 
the wall and we look for a solution of the Navier-Stokes equations which 
satisfies the conditions 

u = w = 0, v = v 0 , y = 0, 

( 2 . 1 ) 


u -*■ 



w w 0 , y , 


where 1 is a length scale in the x direction whilst Uq, Vq and Wq are 
independent velocity scales. We define the parameters A, R, and k by 


A 



R = 


W Q A 


K = V f— ' 1 1/2 

VvU J 


( 2 . 2 ) 


Thus, A is the thickness of the boundary layer associated with the flow in 
the x-y plane. It is known that there exists an exact solution of the 
Navier-Stokes equations corresponding to (2.1). The velocity field of this 
exact solution can be written 


where 


u = u o? u(n) ’ 


v = U of v(n) ’ 


W = W Q w(n ) , 



(2.3) 


(2.4) 


and u, v, w satisfy 

u + = 0, 

v'" + v' 2 - v v" - 1 = 0, 
w"- v w' = 0, 


(2.5) 


v'(0) = 0, v(0) = k , v'(“) = -1, w(0) = 0, w(" ) = 1 
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We note that the solution of the above system has 


v ~ -(n - 6 (k)) when n + », 


where the displacement 6 (k) is, of course, a function of k. We further 
note that if k > 0, then w" will necessarily vanish at some values of n 
so that if u and 7 are, in some sense, negligible in the stability 

calculation we expect inflection point instabilities to occur. 

We restrict our attention to the linear stability of the flow (2.3) to 
disturbances periodic in the z-direction with wavelengths M- . We shall 
take a to be real and look for the corresponding complex frequency 
associated with a for each value of R and k. We perturb the flow (2.3) 
by writing 

“ j[ eu + eRUE] , 


r; — — W + WE, 

“o 

Uq W 

where c - — an d g = exp i a{^- - — ^ — } . It will also be useful to work 

with the quantities cd, X given by 


to - ac, X = -iaRc. 


(2.7) 


The pressure perturbation corresponding to (U, V, W) is pW^ PE and U, V, 
W and P are functions of the nondimensional variable q defined by 
(2.4). The x-dependence of the disturbance which 


we have assumed corresponds 
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to that used by Hammerlin (1955) and enables us to find a solution of the 
linear stability equations by solving ordinary differential equations. Such a 
simplification does not happen in more general boundary layer centrifugal 
instability problems and a self-consistent linear stability analysis leads to 
a partial differential system (see Hall (1983)). It is a routine matter to 
substitute the disturbed flow velocity and pressure fields into the Navier- 
Stokes equations and linearize in U, V, W to obtains 

{_!_ - a 2 }U = -iaRcU + 2u U + Vu' + vlT + iaRw U, (2.8a) 

dn 2 

{ d - « 2 }V = RP' - iaRcV + vV' + v'V + iaRwV, (2.8b) 

dn 2 

2 

( d , - a 2 }w = iaRP - iaRcW + vW' + RVw' + iaRwW, (2.8c) 

1 , 2 1 
dn 

U + V' + iaW = 0. (2.8d) 

The pressure perturbation P and the z component of velocity W can be 
eliminated from (2.8) to give the following pair of coupled equations to 
determine U and V; 

{M + iaRc}u = 2uU + "u'V + "vLT + iaRwU, (2.9a) 


{M + iaRc}MV = iaRwMV - iaRw"V + vMV' + v'MV - 2u'U - 2uU' - u"V - u'V'. 


(2.9b) 
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Here the operator M is defined by 



> 


and we see that (2.9) with R — 0 reduce to the usual equations governing 
the stability of stagnation point flow to Taylor-Gortler instabilities. 
Alternatively by setting u = v = 0 in (2.9b) we obtain the Orr-Sommerfeld 
equation for V. However, it is clear that there is no rational approximation 
to (2.9) which decouples the equations in order to produce the latter 
simplification. Before discussing the numerical solution of (2.9) we must 
derive the appropriate form of the disturbance for n » 1, this will enable 
us to generate boundary conditions at n = n » 1. 

OO 


3. THE BOUNDARY CONDITIONS TO BE IMPOSED AT n = 00 AND THE CONTINUOUS 
SPECTRUM 

Suppose next that we let n 00 in (2.9) and replace u,v,t7 by their 

asymptotic expansions for n » 1 to obtain 

[M - X - 2 - iaR]U + [n - 6 ]U' = 0, (3.1a) 

[M - X - iaR]MV + [n - 6 ]MV' + MV + 2U' = 0, (3.1b) 

where X is given by (2.7). If we look for asymptotic solutions of (3.1) 
with U and V having an algebraic dependence on [ti — 6 ] we obtain the two 
independent solutions 
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2 2 
it f ^-yX+2+iaR-kx „ ✓ ^ x X+l+iaR-Kx 

a. u ~ (ti ) , v ~ (n ) , 


, „ „ / , x X-l+iaR-hx 

b. u = o, v ~ (n ) , 


where n' = n - S. In addition, there exist the following three independent 
exponentially decaying solutions: 


c. U = 0, 


V ~ 


.-an 


> 


d. U = 0, 


, ,v-2-a 2 -X-iaR -n' 2 /2 

(*1 ) e 9 


, , x -3-a 2 -X-iaR -ri' 2 /2 „ [3 + a 2 + X + iaR] U 

e. U ~ (n ) e , V 5 o • 

n' 

2 

Suppose that X is constrained to lie in the half plane X^ < -2 - a 
and we impose the condition U, V, W + 0 when t| + We see that there are 

five independent solutions of (2.9) which have this property and they 
correspond to the asymptotic forms (a) — (e). Each of these solutions can be 
multiplied by an arbitrary constant, a ,**»,a , and the solutions must then be 
combined so as to satisfy 

U(0) = V(0) = V'(0) = 0. (3.2) 


9 

Thus for each X < -2 - a there are two independent eigenfunctions 

r 

corresponding to say, a =0 and a, = 0, We refer to these eigenfunctions 

a d 

as the B and A eigenfunctions respectively. We further note that the 
pressure perturbation corresponding to these eigenfunctions also tends to zero 


when T) + <* 
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Suppose next that X^ is constrained to be in the interval 
2 2 

-2 - a < X r < 1 - a so that (a) no longer gives algebraic decay for U 
when n However, the asymptotic form (b) still corresponds to algebraic 

decay of the disturbance when n - so that there exists one eigenfunction 
(i.e., the B eigenfunction) corrersponding to each X in this strip. In 
particular, if we set X = 0, R = 0, we obtain the continuous spectrum 
0 < a < 1 which was discussed by Hammerlin (1955). 

The eigenfunctions of the continuous spectra discussed above have the 
property that the velocity fields tend to zero algebraically when r) if 

this condition is replaced by one which requires that U, V, W should tend to 
zero exponentially when n + «, then only (c), (d), and (e) can be used to 
generate independent solutions of (2.9). In this case the three independent 
solutions of (2.9) must be combined to satisfy (3.2) and then only a discrete 
spectrum will exist. This possibility was investigated by Wilson and Gladwell 
(1978) who considered the plane stagnation point problem R = 0. The latter 
case is that discussed by Hammerlin (1955) and Wilson and Gladwell found a 
damped discrete spectrum corresponding to disturbances which decay 
exponentially when q + ». It was suggested by the latter authors that the 
continuous spectrum is not physically relevant for the reasons we discuss 
below. 

Suppose that the basic flow which we are considering is the local flow 
near the attachment line on, say, an infinitely long swept cylinder. The 
disturbances which we are calculating must then be matched with a disturbed 
flow in the appropriate 'outer region' unless they decay exponentially when 
11 * in which case the disturbed flow is confined to the boundary layer. 
Thus the physical relevance of the continuous spectrum in this more realistic 
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flow is not known until the disturbance flow field in the 'outer region' has 
been found. It seems that, as yet, the latter structure has not been 

determined so that the relevance of the continuous spectrum remains an open 
question, Wilson and Gladwell argue that the algebraic decay of the 

continuous spectrum eigenfunctions will lead to inconsistencies in the 

matching procedure between the inner and outer region. This is indeed 
possible but perhaps a more likely outcome is that only the unstable part of 
the continuous spectrum would be ruled out by such a procedure. Without 
solving the outer problem the above arguments are merely speculative but it is 
clear that the exponential decay of the discrete spectrum removes any doubt 
about its relevance to more realistic flows. For that reason we shall confine 
our attention to the discrete spectrum but we must bear in mind that the 

continuous spectrum may play an important role the linear stability problem 
under investigation. Indeed, Dhanak and Stuart^) who investigated the 
problem with R = 0, argue that even some of the algebraically growing 

solutions of (2.9) are physically relevant. However, the close agreement 

between experimental observations and our results suggests that the discrete 
modes with exponential decay at infinity dominate the boundary layer stabilty 
characteristics on a swept cylinder. Finally, we note that if U and V are 
to decay exponentially when n + 00 then for sufficiently large values of 
ti, U ~ e ^ ' , V ~ e” 0 ^. Thus the normal velocity component V decays to 

zero more slowly than U. 

Dhanak and J. T. Stuart, 1983, Imperial College, England, personal 
communication, to be published. 
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4. THE NUMERICAL WORK 

The first step in the calculations was to integrate (2.5) by a shooting 
method by guessing an initial value for v"(0) and iterating until 
v'(°°) = -1. The results obtained were in excellent agreement with those 
given in Rosenhead (1963). 

The linear stability equations (2.9) were then integrated subject to (3.2) 
together with the condition that (U, V, W) should tend exponentially to zero 
when h 00 . The first method used was a shooting procedure which though 
able to reproduce Wilson and Gladwell's results at R = 0 without much 
difficulty or expense proved prohibitively expensive at higher values of the 
Reynolds number where instability is possible. For that reason the following 
independent schemes were implemented. 

I. The Riccati Transform Method 

A detailed account of this method can be found in, for example, Aziz 
(1973), Scott (1973), Davey (1977) or Sloan (1977). Davey (1977) has compared 
the efficiency of the method to an orthonormalization method for Orr- 
Sommerfeld equations whilst Wilson and Gladwell (1978) used the method for the 
plane stagnation point problem. The system (2.9) is written in the form 

d Xl 

dT = A lh +A 2 12’ 

d Zz 

*r " a 3 *i + a 4 

where Aj, A 2 , A 3 and A^ are defined by 


(4.1a) 

(4.1b) 
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A ! = 


A 3 " 


and 


0 

-2u‘ 


A 4 " 


1 — 
o 

o 

o 


t— ‘ 

O 

o 

0 0 1 

A 2 = 

0 0 0 

0 0 0 


0 1 0 

L_ _ 

iaRw 

u' 

0 


+ 

i 

1 

a^ X + ia^ Rw" + 

2 — v — , 

a v + u 

V 

0 

0 " 

0 

0 

1 

■2u 

2a ^ + X + v + iaRw 

V 

- (U,V,V') 

|T > y 2 - <u',v~, 

V'"). 


2 - 

-a v - u 


The matrix T, defined by y^ = T jr , satisfies the equation 


£ - A, T - TA 4 + A 2 - TA 3 T, 


(4.2) 


whilst S = T ^ satisfies 


dfi- 55 A 4 S - SAj + A 3 - SA 2 S. (4.3) 

Apart from the definition of y^ and y^ the above formulation follows 
closely that given by Wilson and Gladwell. The no—slip condition requires 
that T satisfies the initial condition 



-15- 


T(0) = 0, (4.4) 

and (4.2) was then integrated from n = 0 to n = rT using a fourth-order 

Runge-Kutta scheme. If p is taken to be large the matrix T contains 

exponentially large-elements so it is convenient to take rT = 1 and calculate 

S before integrating (4.3) from p = p" to n=n »1. 

00 

The three asymptotic forms (c), (d), and (e) found in Section 3 show that 
at large values of rt 



and the eigenvalue X = X(a,R,ic) can be found by iterating to make 
det(SM 1 - M 2 ) = 0. 

The step length of the Runge-Kutta integration and p were varied until 

00 

X was obtained with sufficient accuracy. We now describe an alternative 
method which was used to solve the eigenvalue problem. 
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XX. A Compact Fourth-Order Finite Difference Scheme 

The scheme which we now describe was proposed by Malik, Chuang and 
Hussaini (1982) and is derived by means of the Euler-Maclaurin formula: 


U — 1 

ly 1 '* _ ^ 


- (^1. + 
2 '■dn 


df 


k-1 


dn 


) 12 i 


^ 2 k 

kfd z r 


d 2 f k - ! 


dn‘ 


dn‘ 


■) + 00O, (4.5) 


where 

'F k ='l , (n k ) and h k = ’Ifc " n k-r 


The nodes are distributed to resolve singular layers so that 


n k = L(k-l)/(Ng - (k-1)), k= 1 ,2 ,• • • ,N + 1, (4.6) 

where N + 1 is the total number of nodes, g = (n^ + the location 
of the boundary layer edge, and L a scaling parameter. We chose L to be 
twice the height (from the wall) of the nodal point at which the 
nondimensional mean velocity w assumed the value 0.5. It was found that for 
a given number of points, this choice of L in the grid distribution function 
(4.6), yielded maximum accuracy. 

In order to apply the above compact difference scheme to (2.9a,b), it is 
necessary to write them as a set of first-order differential equations. They 
are rewritten as 


d<t> ± 

dn 


6 

I 

j=l 


*ij 9 j 


i = 1,2,» • • ,6, 


(4.7) 


where 
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♦l 


= v. 


^2 


d<f>j 

dir * 


*3 




9 


*5 


= u. 



9 


with the boundary conditions 


+ 1 = *2 = +5 = ° 


at n = 0, (A. 8) 


$2 + a<l> i = $3 + ai >2 = <(> 6 + n<f> 5 = 0, n = T) m . (4.9) 

The latter conditions follow from the discussion of Section 3 where it was 
shown that for sufficiently large values of ri 


U 


-n 2 /2 


.-an 


The nonzero elements of the matrix ( a fj) are 


a I2 - *• 


a 23 1( 


a 34 ’ *• 


4 2 3 2 _ 

= -a + ia u)R - ia Rw - a v - iaR w" - u' 


9 
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a 


42 


2 - 
-a v 


iT 


9 


a 


43 



ia)R + iaRw + v', 


a 


44 


= v, 


a 45 ' ■ 2u ' 1 


a 46 ’ " 2u > 


a, c = a - ia)R + 2u + iaRw, 

03 


We set 


a 66 ’ v ‘ 


*- { v- £- 


c > 

J 1 

1 d 2 * _ J 

L 6 . ) 

1& a « **j 

t* d^ 2 i 

(/.i b « 


where 


da , , 6 

b ij " dn + a U a *j 


and thus ( 4 • 5) becomes 


2 

k \ r k k h k y , k k 

*i ' T ^ "i J ♦ j + 12 jij b ij ♦ j 


ii h i 6 

^ 1 + TT- I 


h 2 fi 

k-1 .k-1 . "k r , k-1 k-1 

*J + 17 J i 1 b lJ 


J 


= 0, 


(4.10) 


k = 2,3,» • • ,N + 1. 
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It is possible to write the above equation system along with the boundary 
conditions (4.8) - (4.9) in block-tridiagonal form: 


\ <i> k 1 + B k <j> k + <f, k+1 = H, k = 1 , 2, • • • ,N + 1 (4.11) 


where A k , B k , C k are 6x6 matrices defined below and H is a 6x1 null 

matrix. 


\ 


6 columns 

Nonzero elementsl 3 rows 
0 I 3 rows 


; k = 2,« • •• ,N + 1 , (4.12) 


6 columns 


Nonzero elements 


3 rows 


3 rows 


; k =!,•••• ,N. (4.13) 


B k are full matrices with k = 1 and k = N + 1 having the structure 


B 


1 


6 columns 


E 


Nonzero elements 


3 rows 
3 rows 


(4.14) 


- 20 - 


B 


N+l 


6 columns 
Nonzero elements 
F 


3 rows 
3 rows 


(4.15) 


where E and F are 3x6 matrices representing the bottom and top boundary 
conditions (4.8) and (4.9) respectively. Thus we obtain 


E 


1 

0 

0 


0 

1 

0 


0 

0 

0 


0 

0 

0 


0 

0 

1 


0 

0 

0 


(4.16) 


and 


a 


F = 


0 

0 


1 

a 

0 


0 

1 

0 


0 0 

0 0 

0 n 


0 

0 

1 


(4.17) 


We can write (4.11) in the form 

£♦ = H (4.18) 

where £ = [A k , B k , C k ] . 

Assuming an estimate of the eigenvalue is available, we solve (4.18) 

directly. In order to avoid the trivial solution, nonhomogeneous boundary 

conditions are imposed at the wall. Specifically, the boundary condition 

<)> (0) = 0 is replaced by <(» _ (0) » 1. This is equivalent to normalizing the 
1 J 2 

eigensolution by the value of at the wall. Matrix E of (4.16) now 

dn 


becomes 



- 21 - 


and 


E = 


0 

0 

0 


0 

1 

0 


1 

0 

0 


0 0 
0 0 
0 1 


0 

0 

0 


H = (1, 0, 0,*»» ,0) T . 


The system (4.18) is nonhoraogeneous and the nontrivial solution is 
obtained using block LU factorization. Newton's method is then used to 
iterate on the eigenvalue such that the remaining boundary condition 
<f>j(0) = 0 is satisfied. 

In order to generate neutral curve (w^ = 0), a solution $ is first 

obtained for assumed values of a and w • The corrections Aa and Aw 

r r 

are then determined from the equations: 

Hi ( 0 ) Hi ( 0 ) 

+ 1 (0) Aa + 9 ^ A( 0 r = 0, (4.19a) 

r r 

9<}> 1 ( 0 ) 3 *j ( 0 ) 

<}>, (0) +- 5 — - Aa + ■* — - Aw = 0. (4.19b) 

T l. 9a 3o) r 

i r 


Here <t>j(0) is known from the solution <|> just obtained whilst the 
derivatives with respect to a and w^, are obtained by solving 


and 



P ^4 _ 3eC . 

^ 3w “ 3w 9 
r r 


(4.20) 


(4.21) 
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The process is repeated until <f>^(0) vanishes within preassigned 

tolerances. We see that (4.18), (4.20) and (4.21) can be solved with the same 
LU factorizations and both the eigenvalue and eigenfunction are obtained. 

The two schemes described above were first used to calculate the 
eigenvalues of (2.9) with R = 0. The results obtained were found to be 

consistent with those given by Wilson and Gladwell (1978). At higher values 
of R, where instability is possible, the Riccati method gave consistent 
results provided that enough grid points were used. For Reynolds numbers less 
than ~ 1000 it was found necessary to use at least five hundred steps (with 
r\ m = 10) to obtain a neutral eigenvalue correct to four significant 

figures. In Table I we have illustrated the convergence of the finite 

difference scheme in this region. 

The convergence of the latter scheme is significantly faster than that of 
the Riccati method and the neutral values of a , c^ correct to four 
significant figures can be obtained with only forty-one grid points. 
Moreover, it was found that the iteration to an eigenvalue with the Riccati 
method was much more sensitive to the initial guess for the eigenvalue than 
was the case with the finite difference scheme. In view of the fact that the 
eigenfunctions cannot be calculated directly using the Riccati method because 
of numerical instabilities it was therefore decided to use the finite 

difference scheme to generate the neutral curves and use the Riccati method to 
check selected eigenvalues. 

In Figure 1 we have shown the neutral wavenumber and frequency to = ctc^ 
as functions of the Reynolds number R for k = -.8, -.4, 0, .1. The 
critical values corresponding to the zero suction case < = 0 are 
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R = 583.1, a = .288, etc = .111. 

r 


It xs clear from Figure 1 that suction and blowing have stabilizing and 
destabilizing effects on the flow respectively. The shape of the neutral 
curves for k < 0 is typical of viscous instabilities of flows without 
inflection points whereas the open neutral curves for k > 0 are to be 
expected since w has an inflection point in this case. It should be noted 

that, although (2.9) is not the Orr-Sommerfeld equation, the inviscid limit of 
(2.9) is the Rayleigh equation so that in the absence of an inflection point 
the flow is stable at infinite Reynolds numbers. 

In Figure 2 we have shown some of the eigenfunctions corresponding to the 

neutral curves of Figure 1. The eigenfunctions have been normalized in each 

case such that the maximum magnitude of each velocity component is unity. It 

can be seen that the x-velocity component goes to zero more quickly than the 

other two components. This is, of course, to be expected since we have shown 

2 , 

earlier that for large n, U ~ e whilst the other velocity components go 

to zero like e ari . We also notice that for the cases with k > 0 the 

eigenfunctions have their critical layers further away from the wall than is 
the case with k < 0. 

In Figure 3 have have plotted the critical Reynolds numbers corresponding 

V n 

to k - +.8, +.4, 0, -.1, -.15, -.2. We note that jl- = so that the 

W ° R 

significant rise in the critical Reynolds number corresponding to k = — .1 
, V 0 -4 

corresponds to 10 which is well within the range of practical 

0 

feasibility. Figure 3 suggests that the critical Reynolds number can be made 
arbitrarily large by taking -k << 1. In order to determine whether this is 
the case we shall in the next section consider the structure of the neutral 
curve in the limit k 
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5. THE LARGE SUCTION LIMIT k -► -» 

We first consider the limiting form of the solution of (2.5) when 
K 4- For n ~ 0(1), it is easily shown that the asymptotic form of v is 

2 

v~ k - ^— + ••• , (5.1a) 

2k 

whilst 

u~— +••*, (5.1b) 

k * 

It follows from the z momentum equation that w then approaches unity on 
the length scale |k| - 1 so that if we write 

5 = n |k [ 


then for £ ~ 0(1) 

u ~ — y v ~ k - „ + w= (1 — )+••*. (5.2) 

K^ 2k J 

Thus w is given to first order by the asymptotic suction velocity profile 
and before determining the corresponding asymptotic structure of U, V we 
first write 

a = |k | a, (5.3a) 

R = |k|r. (5.3b) 

The above rescaling is suggested by the fact that the appropriate length 
scaled for the instability is now v/V Q and not A. 


It now remains for us 



-25- 


to rewrite (2.9) in terms of £ and expand U,V in the form 


U 0 U i 

u= 3 + 4 + * 

K K 


V = V n + — + 

0 K 


where U Q , Vj, etc. are functions only of % . It is then a routine matter to 
substitute the above expansions into (2.9) and determine U Q , V 0 , etc. We 
find that Vq satisfied the equation 


|^2 ~ a 2 + iaRcj jij-g. -a 2 jv Q « iaR[i-e" ? ]v () + iaRe < V 




0* 


(5.4) 


which is just the modified Orr-Somraerf eld equation for the asymptotic suction 
profile. The appropriate boundary conditions for (5.4) are 


v 0 = V o = °> 5=0, V Q > 0, 5 + ». (5.5) 

There is no continuous spectrum of (5.4) corresponding to (5.5) and the 
critical Reynolds number associated with the discrete spectrum is given by 
Hocking (1974) as 

R = 54,370 

whilst the critical values of a, c are 
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a = .1555, c = .150. 

Hence, we see that in the limit of large suction (2.9) simplifies to the 
modified Orr-Somraerfeld equation for the asymptotic suction profile and the 
critical value of R is given asymptotically by 

R ~ 54,370 |ie|. (5.6) 


The massive stabilizing effect of suction suggested by Figure 3 is therefore 
not surprising. It follows from (5.6) that R can be made arbitrarily large 


by taking the limit tc + -» but that instability then always occurs when 

— 2 . = - ■ and is confined to a layer of thickness (v/V n ) at the wall. 

Wq 54,370 u 


6. DISCUSSION 

We have shown that the attachment line boundary layer flow (2.3) is 
susceptible to travelling wave instabilities which propagate along the 
attachment line. These disturbances decay exponentially at the edge of the 
boundary layer and at zero Reynolds number correspond to the eigenvalues of 
the plane stagnation point stability problem found by Wilson and Gladwell 
(1979). The exponential decay of the eigenfunctions ensures that in more 
realistic situations there will be no difficulty matching the solutions onto 
the appropriate outer flow field. We now turn to the relevance of our work to 
the available experimental results which correspond to the flows around swept 
flat— nosed wings and circular cylinders. Before doing so, it is of course 

necessry to indicate the relevance of our basic flow to these 
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configurations. A discussion of the boundary layer flow on an infinitely long 
swept cylinder is given in Chapter VII of Rosenhead (1963) and is based on the 
work of Schubart (1945, unpublished), and Sears (1948, 1952). 

The procedure is essentially identical to that appropriate to the two- 
dimensional case and outside the boundary layer the basic flow is expanded in 
the form 


In the boundary layer the velocity components are expanded in powers of — 
and the first order flow is exactly that which we have described in Section 2. 
Thus, sufficiently close to the attachment line, the basic flow which we have 
calculated in Section 2 gives a good approximation to the attachment line 
boundary layer on a swept cylinder. The range of validity of our theory will 

therefore depend on the cross-sectional shape of the cylinder and the Reynolds 
number. 

The experimental investigations most relevant to the present work are due 
to Gaster (1967), Pfenninger and Bacon (1969) and Poll (1979, 1980). The 
first two investigations were concerned with the attachment lines on swept 
wings whilst Poll performed experiments on a swept circular cylinder. 
However, the wing section used by Pfenninger and Bacon had a flat nose so it 
is likely that our theoretical work is most relevant to their experiments. It 
is well-known that the attachment line is turbulent when large amplitude 
disturbances are present. These large amplitude disturbances convected into 
the attachment line can be prevented by various means so that the linear 
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critical Reynolds number almost certainly gives an upper limit to the laminar 
flow regime. Gaster found that the attachment line boundary layer is stable 
up to R ~ 420 which was the maximum value of the Reynolds number which could 
be achieved experimentally. 

Pfenninger and Bacon (1969) and Poll (1979, 1980) measured equilibrium 
frequencies of naturally occuring small disturbances at various values of the 
Reynolds number. We have indicated these experimental results in Figure 4 
where we have also shown the neutral curves without suction. The experimental 
points lie close to the lower branch of the neutral curve and the critical 
Reynolds number is close to the minimum Reynolds number at which Pfenninger 
and Bacon measured small amplitude equilibrium disturbances. It can be seen 
in Figure 7 of Pfenninger and Bacon (1969) that by introducing large amplitude 
disturbances into the boundary layer by means of trip wires it is possible to 
induce equilibrium disturbances at significantly lower values of the Reynolds 
number. The latter authors also give experimental points in the presence of 
suction but it is not clear from the results given whether the flow was 
stabilized or destabilized by the suction. The amount of suction used 
experimentally is not given by Pfenninger and Bacon but, in view of the 
results shown in Figure 3 of this paper, we assume that it was extremely 
small. 

It seems likely that the instability mechanism which we have discussed is 
responsible for the disturbances measured experimetnally • Some further 
evidence for the relevance of our model is provided by Figure 5 which is taken 
from Poll (1980). Poll used two hot wires at the same spanwise location but 
different chordwise locations to measure the time history of the spanwise 
disturbance velocity component. The figure suggests that the instability is 
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indeed two-dimensional and takes the form of a time-ha rraonic disturbance 
propagating along the attachment line. 

Finally we close by making a few comments about the extension of our 
analysis to include the effect of nonlinearity. Surprisingly, the disturbance 
structure (2.6) can still be retained and enables us to take out the x- 
dependence of the problem even in the nonlinear regime. Thus it is possible 
in principle to calculate the weakly nonlinear development of the disturbances 
considered in this paper. Such a calculation might explain the origin of the 
subcritical equilibrium disturbances measured by Pfenninger and Bacon (1969). 
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Table I. The neutral eigenvalues as a function of the number of grid 


points of 

the finite difference scheme for 

R = 800, k = 0. 

10 

.3300581 

.1226919 

20 

.3378719 

.1267951 

40 

.3384238 

.1270776 

80 

.3384613 

.1270965 

160 

.3384638 

.1270977 
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Figure Captions 


Figure la : 

The neutral curves in the a - R plane for different values of k. 

Figure lb : 

The neutral curves in the u> - R plane for different values of k. 
Figure 2: 

The neutral eigenfunctions for: (a) R = 583.4, a = 28, 

(b) R = 4000, a = .286, (c) R = 44., a = .4014, (d) R = 119, 

a = .3587, (e) R = 979, a = .285. 

Figure 3: 

The dependence of the critical Reynolds number on the suction parameters. 
Figure 4: 

A comparison between theory and experiment for k = 0. 

Figure 5: 

The spanwise disturbance velocity components measured by Poll (1980) at 
neighbouring values of x with z fixed. 






Figure lb 
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Figure 3 
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Figure 4 


O Pfenninger and Bacon (1969) 
□ Poll (1979, 1980) 
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